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Рассматривается и решается начально-краевая задача по математическому моделированию поля концен- 
трации электролита в управляемом электрохимическом сопротивлении (УЭХС) с плоскими электродами. 
Модель позволяет вести полное исследование электрохимических процессов во всём объёме УЭХС. Выпол- 
нены контрольные расчёты пространственно-временного распределения концентрации при постоянном 
управляющем токе УЭХС. 

Ключевые слова: краевые условия, параболическое уравнение диффузии, плотность электрического 
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Введение. УЭХС представляет собой цилиндрический проводник длиною / прямоугольного се- 
чения (подложка) ахг , изготовленный из стеклографита. Если на боковую поверхность УЭХС 
наносить электрохимическим способом тонкий слой металла, то сечение слоёного проводника 
будет увеличиваться и его проводимость возрастёт. Распределение слоя осаждаемого металла 
по боковой поверхности УЭХС полностью определяется распределением концентрационного 
поля электролита. 

В статье рассматривается постановка и решение начально-краевой задачи для про- 
странственно-временного распределения концентрации электролита в УЭХС. 

Основные допущения: 

— лимитирующей стадией электродных процессов в цепи управления УЭХС является мо- 
лекулярная диффузия; 

— коэффициент диффузии электролита есть величина постоянная и не зависит от тем- 
пературы; 

— плотность электрического тока по поверхности управляющего электрода распределена 
равномерно; 

— пространственно-временное распределение концентрации электролита является плос- 
ко параллельным и зависит от координат х; уи времени Е 

— толщина наносимого металла значительно меньше поперечных размеров УЭХС. 
Постановка задачи. Согласно принятым выше основным допущениям, все физические поля в 


УЭХС являются плоскопараллельными и не зависят от координаты 2 (рис. 1), перпендикулярной 
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плоскости хОу. Слой металла, наносимый на подложку УЭХС, является ничтожно тонким по 
сравнению с геометрическими размерами самой подложки и толщиной слоя электролита. 
На основании 2-го закона Фика относительно концентрационного поля электролита 


С (х; у; В ставится следующая начально-краевая задача. 








п ртоя, еоуч ” 
Се и:0) =. (2) 

оС кои) - М, в) 
в) М (жуй), (4) 

2 (6;;2) =0, (5) 

(ви) = 0, (6) 


где р — коэффициент диффузии электролита; С — начальная концентрация электролита; а — 


толщина УЭХС; 5(х;Е) — неизвестное и подлежащее определению распределение плотно- 


сти тока на подложке; /М— электродно-кинетическая константа. 





Рис. 1. Схема включения УЭХС 


Задачу (1) — (6) удобнее всего решать операторным методом Лапласа. Для этого введём 


следующие соответствия: 
С(жу:#) = С(х;у:рР)} БЬЮ =Ь(Р); 5(х;у;Ё) .= 9(х;у;р). 


В пространстве изображений относительно С (х;у;р) имеем краевую задачу: 





(7) 


се те 
9х? ду?) 


рС(х;у:Р)- С. -2 ы 
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2 кур) М), () 
2 (ину) = Мур), (9) 
2 (зу;р) 0, (0 
2 (ррузр) 0. ив 


Преобразуем уравнение (7) к виду: 





2/- 2/- , 
ИЯ РС Со. 
0х’ ду’ О р 





(12) 


Краевую задачу (8) — (12) будем решать методом разделения переменных. Для этого 


введём новую функцию и (х;у;р) по правилу: 








. : С 
и(х;у;рВ) о (13) 
Для неё получается вспомогательная краевая задача: 
сы: р. 
9х? г ду? В“ 0, (14) 
20 озр) - МР) (15) 
бу Г а/ ' 
су(иур)=М8(ж:р), (16) 
ину 
5х (6У:Р)=0, (17) 
ди 
(/:и: =0. 
| ;У;Р)=0 (18) 


Общее решение уравнения (14) определяем как суперпозицию частных режимов [1]: 


и (х;у;р) -[м (р) Бу +М бодек ут 
+ >", (р) 5+ му +М, (р)сН + ду. (х), 


где р, = и, К=1,2,... — собственные числа краевой задачи (14) — (18), 


(19) 








—> 





Х,(х)= есь т», К =1,2,... — собственные функции, М, (р), № (р), М, (р), №, (р) — 


произвольные постоянные. 
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Для конкретного определения коэффициентов разложения в (19) продифференцируем 


2убву:Р)- м о (р рев РУз (р зу 
"ним „рае ом рав Ву 69 


Из (20) при у= биу= РВ соответственно получаем: 


| В 
ине) Мрт чм, (ох, © 


ифитр) = [В мови (р) РЯ т 
+2. Би [доза ещь м, (65 ня же 


Подставим теперь (15) в (21), а (16) в (22): 


ма (2) _ ры м, (в) ям, (р), (23) 
мокоуаь Ве ньдрунь т 
+У ть [м м (реп +щи+м, фоуан]х, (х). 


К= 
Применяя к (23) стандартную процедуру поиска коэффициентов разложения Фурье, по- 


его по у: 


(20) 


(22) 


$ 





Мб (х;р)= 


я < 


(24) 


: 


лучаем: 
М, (р) 
М (р) т [2 (25) 
а\! - | 
М, (2)=0,К=1,2, ... (26) 


Разложение (24) с учётом (25) и (26) принимает вид 


Гу. Р, 1х Р р 
№ (х;р) у 2-м м сям, рев |2 о ГМ (р)58 [чи ИХ, (х). (27) 
Из (27) с использованием выше приведённой процедуры поиска коэффициентов имеем: 


В, (Рр)- аа 
РЕ 
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М Ви (р 
к Е 
где неизвестные коэффициенты разложения В, ( р), В (р) в ряде Фурье неизвестной плотности 


тока 5(х;р) определяются соотношениями: 


= [6ббр)тх, в, (2)- [869р)х, (р)ах. 


Эти коэффициенты могут быть найдены только в результате совместного рассмотрения 
концентрационного и электрического полей в УЭХС, что будет сделано далее. 


Подставим найденные коэффициенты М, (р), № (р) и М, (р) в выражение (19). 


1 
о 


и (х;у; р) = ре =. Е ав) У) ы 
.Я МВ, (р) а Рщу-х, 


и. ты аи |+ 











Подставляя выражение (28) в (13), находим аналитическое описание концентрационного 


поля электролита. 


@ т 


ЩЕ 
"2 тя а Р+щу.Х, (*) 
Ре Руль 


С помощью полученной формулы (29) можно на данном этапе получить явное выраже- 








С(х;у;р)= 








ние для расчёта концентрационного поля лишь при равномерном распределении плотности то- 
ка на подложке, когда 


Г, (р) 


==) (30) 


5(х;р) = 





а разложение этой плотности в ряд Фурье 


1 Р) - в, (р), 


приводит к таким результатам: 


В, (р) ==—^, (31) 


Вестник ДГТУ. 2012. № 1 (62), вып. 2 








В (2) =0, Кд (32) 
Подстановка (31) и (32) в (29) окончательно определяет концентрационное поле элек- 


тролита при равномерном распределении плотности тока: 


а М В 
вЫ ыы > =. 


аа 
р р 2 





При включении входной цепи УЭХС на постоянный ток Г, 


1 (2). (34) 


Подставим (34) в (33). 
С 


: М р р 
Су; р) =— +1 5И } 
(у;Р) т - р Я (35) 
ар Рен. |Р 7 
р р2 


Для обращения С (у;р) рассмотрим вспомогательные функции: 


На 





Е (у: р) = ——^, (36) 
Ра ее 
р р2 
| Е (у; 
Е, (у;р) = и й (37) 


Для поиска оригинала # (у;#) .= Е (у; р) найдём вначале полюса А, (у;р). Очевидно, 


что число р = 0 не является таким полюсом. Эти полюса определяют из условия 


|2 о (38) 
Б2 


Для решения уравнения (38) введём новую переменную 27 по правилу 


|= й - №, (39) 


где ] = /-1. Подставим (39) в (38). Получим уравнение сп} =0, которое равносильно урав- 
нению со52 =0. 


Корни последнего уравнения имеют вид: 


2, = (2+1), ВЕ. (40) 


Подстановка (40) в (39) определяет полюса 2. (у;р): 


(2к+1) РР 
р? 








р, = ‚ КЕЙ. (41) 
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Далее находим вычеты 


а Ве | 
. о 5] ‚о тат. о в а 


р п В п В 
ее : (2К ту - "| я ыи[ (2% + АЙ Е "] окно, 


— = е 1? 


: ь 5" (2к +107) 


вк) 


Оригинал # (у;Ё) находят по схеме 




















(уз) - вез {е”А, (у;р)} (42) 


Запись (42) в развёрнутом виде даёт такой результат: 


: р 
ни + Эту = "] екяуто, 
е 


40 - 2 1 (43) 
0 ; Ш 
т (2^ +1)> 


В 





[М 


ин 


> 
И 


Оригинал Ь, (у;Ё) .= Е (у;р) находят, интегрируя (43) по времени в пределах от 0 до & 


. ни (6 + Эту _ "| ака то, 
в ре | (44) 


г? (2к +1) яп(2к +1) 





з|> 


(уй)= 


Оригинал концентрации С (у;Ё) связан с № (у;Ё) соотношением 


С(у;#) =С +1 Ву У;Ё). (45) 


Подставим (44) в (45). 


в 
р ни (2+1) [у- | окне, 
С(у;Ё)=С и. И 1 е " (46) 


ы (2К+1) 5 п(2к +1)5 








Полученное выражение для С(у; #) позволяет с любой степенью точности построить по- 
верхность изменения концентрации электролита в УЭХС. 

При расчётах по последней формуле бесконечный ряд заменяется конечной суммой. 
Число слагаемых М в ней задаётся погрешностью расчёта концентрации А. Обозначая разницу 


между точным значением концентрации и приближённым СИМВОЛОМ Ки ‚ получим: 
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. п р 
АМ В ни (26 с я = "| (кт, 
ети 1-е " | (47) 
ИЯ Ок п (2 +1) 
Потребуем, чтобы 
|| < А. (48) 
Рассмотрим вспомогательный ряд 
4ВМ _ < 
5102. 
пга/ 2 (2к т 
ы - и К„ существует соотношение [К„|<О,. 





Между его остатком О„ = 5 


п?а/ 2 (2К+ р. 


Поэтому, если выполнить условие О, < Д, то с гарантией будет выполнено и требование (48). 


Для остатка О„ справедлива интегральная оценка [2]: 


4ВМ _ 4ВМ 1 
1 : 
©» < па/ [5 у5 т пга/ °М 





Число М определяется как решение неравенства 
р Г. г <А. 
При исходных данных 1 = 0,01 м; а = 0,015 м; / = 0,08 м; 4 = 0,2 А; С, =3 кмоль/МЗ, 
№М= 5,45 (кмоль/м3)м"А-", Л =1.103 кмоль/м? были рассчитаны концентрационные профили 


С (у 6) (см. рис. 2) и временные зависимости С (у; 8 (рис. 3) соответственно при фиксирован- 


ных ви и. 





м 


ы 0, 0025 0,005 0, 0075 


Рис. 2. Концентрационные профили электролита 


23 


Физико-математические науки 











рэ 50 29 


Рис. 3. Временные зависимости концентрации электролита 


Заключение. Полученная математическая модель концентрационного поля электролита в 
УЭХС успешно может быть использована для расчёта режимов электрохимического осаждения 
металла на его боковую поверхность. Сопряжение этой модели с описанием электрического по- 
ля позволяет получить математическую модель УЭХС как элемента электрической цепи. 
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МАТНЕМАТТСАЕ 5ТМОТАТТОМ ОЕ ЕТЕСТКОГУТЕ СОМСЕМТКАТТОМ ЕТЕГО ТМ 
СОМТВОГЕАВЕЕ ЕТЕСТКОСНЕМТСАЕ ВЕЗТЗТАМСЕ 


Т. Р. ЗКаКипома, У. У. бегазйтепко 
(Роп Зае Теспткса! Упмег®Ку) 


Тре та! Боипаагу угие ргоМет оп 1е та етайса/ ути!айоп оЁ {те еесвоу сопсепгаНоп Пе т те 
сотгойаЫе еесосвет/са/ гея ${апсе (СЕК) ийН риапе еесосдЕс 15 сопуаегеа апа зо/№е4. Тре тоде! реттй5 
{0 сопаисЕ а сотр гезеагсй оЁ 1е еесосветгса! ргосез5ез т Ви-Б/азЁ СЕК. Тре свеск сотриаНоп оЁ 
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